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Abstract. We study the present, flat isotropic universe in modified gravity. 
We use the Palatini (metric-affine) variational principle and the Einstein (metric- 
compatible connected) conformal frame. We show that the energy density scaling 
deviates from the usual scaling for nonrelativistic matter, and the largest deviation 
occurs in the present epoch. We find that the current deceleration parameter derived 
from the apparent matter density parameter is consistent with observations. There 
is also a small overlap between the predicted and observed values for the redshift 
derivative of the deceleration parameter. The predicted redshift of the deceleration- 
to-acceleration transition agrees with that in the ACDM model but it is larger than 
the value estimated from SNIa observations. 
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1. Introduction 

The most accepted explanation of current cosmic acceleration jU El EI is that the universe 
is dominated by dark energy However, it is also possible to modify Einstein's general 
relativity to obtain gravitational field equations that allow accelerated expansion. A 
particular class of alternative theories of gravity that has recently attracted a lot of 
interest is that of the f{R) gravity models, in which the gravitational Lagrangian is 
a function of the curvature scalar R 5 . It has been shown that cosmic acceleration 
may originate from the addition of a term R~^ (or other negative powers of R) to the 
Einstein-Hilbert Lagrangian P EJ CHI HI] . 

As in general relativity, these models obtain the field equations by varying the total 
action for both the field and matter. In this paper, we use the metric-affine (Palatini) 
variational principle, according to which the metric and connection are considered as 
geometrically independent quantities, and the action is varied with respect to both of 
them 1121. The other one is the metric (Einstein-Hilbert) variational principle, according 
to which the action is varied with respect to the metric tensor g^^, and the affine 
connection coefficients are the Christoffel symbols of g^,y. Both approaches give the 
same result only if we use the standard Einstein-Hilbert action. The field equations in 
the Palatini formalism are second-order differential equations, while for metric theories 
they are fourth-order ^Hj- Another remarkable property of the metric-affine approach is 
that the field equations in vacuum reduce to the standard Einstein equations of general 
relativity with a cosmological constant |14j . 

One can show that any of these theories of gravitation is conformally equivalent 
to the Einstein theory of the gravitational field interacting with additional matter 
fields jin]. This can be done by means of a Legendre transformation, which in 
classical mechanics replaces the Lagrangian of a mechanical system with the Helmholtz 
Lagrangian fH]- For an f{R) gravity, the scalar degree of freedom due to the occurrence 
of nonlinear second-order terms in the Lagrangian is transformed into an auxiliary scalar 
field 0. The set of variables {g^u, 4>) is commonly called the Jordan conformal frame. 

In the Jordan frame, the connection is not metric compatible. The compatibility can 
be restored by a certain conformal transformation of the metric: g^i, h^j^^ = f'{R)g^^. 
The new set {hfj,^, 0) is called the Einstein conformal frame. Although both frames are 
equivalent mathematically, they are not equivalent physically J7] , and the interpretation 
of cosmological observations can drastically change depending on the adopted frame |18j . 
In the Palatini formalism, which frame is physical is a matter of choice, although this 
question should be ultimately answered by experiment or observation. 

The consistence of metric-affine f{R) gravities with cosmological observations has 
already been studied for the Jordan frame ^Hl- In this paper, we regard the Einstein 
conformal frame as the physical one. We assume that the gravitational Lagrangian 
contains the usual Einstein-Hilbert part R and the 1/R term which causes current 
cosmic acceleration [7 . In section 2, we review the metric-affine formalism for an f{R) 
gravity in the Einstein frame [201 • section 3, we apply the gravitational field equations 
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for the R + l/R Lagrangian to a homogeneous and isotropic universe, and solve them for 
the early and late universe in the quadratic approximation. In section 4, we derive the 
relations between the present deceleration parameter qo, the matter density parameter 
Qm and the redshift of the deceleration-to-acceleration transition Zt in the Einstein 
frame of the metric-affine l/i?-modified gravity. We also compare these relations with 
cosmological observations. The results are summarized in section 5. 

2. Metric-afRne formalism in the Einstein conformal frame 

The action for an f{R) gravity in the original (Jordan) frame with the metric g^^ is 
given by 

Sj = -T^ I d'x[^gL{R)]+Sm{9,u,i^). (1) 

Here, ^/—gL{R) is a Lagrangian density that depends on the curvature scalar R = 
R^yiV p^)g^^ , Sm is the action for matter represented symbolically by ip and independent 
of the connection, and k = Tildes indicate quantities calculated in the Jordan 

frame. 

Variation of the action Sj with respect to (jf^i, yields the field equations 

L'{R)R^, - ^L{R)g^, = kT^,, (2) 

where the dynamical energy-momentum tensor of matter is generated by the Jordan 
metric tensor: 

5Sm = ^ I d^xJ^gT^jr". (3) 
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and the prime denotes the derivative of a function with respect to its variable. From 
variation of Sj with respect to the connection F^^, it follows that this connection is 
given by the Christoffel symbols of the conformally transformed metric pUj : 

g^,^ = L'{R)~g^^. (4) 

The metric g^^^ defines the Einstein frame, in which the connection is metric compatible. 
The action (0) is dynamically equivalent to the following Helmholtz action [T^ ll6j: 

1 ^.4. 



= P- (6) 

R=4,(p) 



Sh = -— d^x^~~g[L{<p{p))+p{R-<p{p))] + SU~9,.,^), (5) 
where p is a new scalar variable. The function (f){p) is determined by 
dL{R) 
OR 

From equations (jH) and (jB)), it follows that 

= i?L'(0). (7) 

In the Einstein frame, the action (0) becomes the standard Einstein-Hilbert action 
of general relativity with additional scalar field: 

S. = -^ /AV^[fl-4^ + «l+S„(p-V,^). (8) 



2kc J p p'^ 
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Here, R is the curvature scalar of the metric g^^. Choosing as the scalar variable leads 
to 120] 

Se = -i^J rf^xv/^[i?-2r(0)] + S'™([L'(0)]-V,^), (9) 
where V^(0) is the effective potential, 

<PL'{4>) - m 



Variation of the action with respect to g^^, yields the equation of the gravitational 
field in the Einstein frame ^Ul 1^ '■ 
1 kT 

Rt^u - ^Rg,^u = - V{(f))g^^, (11) 

while variation with respect to reproduces equation (|7j). Contracting (fTTj) with the 
metric tensor g'^'^ gives 

kT 

which is an equation for R since both and T = T^ud^'^ depend only on R due to ^ 
and ((121). This is equivalent to [2(3j 

0L'(0) - 2L(0) = KTL'{(j)). (13) 
Let us consider the CDTT Lagrangian [3 |HI 

^(0) = 0-^, (14) 
where a is a positive constant. t| Equation (fT^ for T = yields one de Sitter solution, 

4>ca = -tt, (15) 

describing present cosmic acceleration [S]. The corresponding value of the curvature 
scalar in the Einstein frame is determined by equation (|7|): 



^ = TT^' (16) 
from which we obtain 

Rca = (17) 

Finally, the Einstein equations (fTT|) become 

1 i^Tf,^ S 

R^lu - T^Rg^u = ofi 77~1^^3^lu■ (18) 

Equations (fT^ and (fTH)) give the relation between matter and geometry in the Einstein 
frame of l/-R-modified gravity. § 

t Equation lO states that (f) is the curvature scalar in the Jordan frame R. The Lagrangian p4|l 
can thus be written as L{R) = R ^ sjid such a model is referred to as the R + 1/R gravity or 
l/i?-modified gravity. 

I The Lagrangian of has an additional term R'^ , which is the simplest way of introducing inflation 
in a metric-affine f{R) gravity [221 . In this work, we omit this term since we are interested in later 
epochs. 

§ The value of a is on the order of 10^^'^m^^ 
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3. The field equations in the early and late universe 

We begin with the Einstein equations for the early universe in which the 1/R term does 
not dominate. In such a universe kT ^ a and we can study the field equations in the 
second approximation of a small quantity Equation (|T^ becomes cubic in 0: 

03 + kT(I)^ - a'cj) + = 0, (19) 

and its solution which deviates from —nT (the solution for a = 0) by terms linear and 
quadratic in is 

4>=-KT-^ + 0{a'). (20) 
In this approximation, the Einstein equations become 

R^iu — ^RQhu = l^effT^u + ^effOfiu, (21) 
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where 

^eff = k(i - (22) 
is the effective gravitational coupling constant and 

A.// =5^ (23) 

is the effective cosmological constant. 

We now proceed to the field equations in the early Friedmann-Lemaitre-Robertson- 
Walker (FLRW) universe. Let us consider a homogeneous and isotropic universe which 
is spatially fiat In this case, the interval is given by 

ds^ = c^dt^ - a^{t){dx^ + dy"^ + dz^), (24) 

where a{t) is the scale factor. Moreover, the energy-momentum tensor of matter in the 
comoving frame of reference is diagonal: 

T; = dtag{e, -P, -P, -P), (25) 

where e is the energy density and P denotes pressure. For this metric and for the case 
of dust (P = 0), equation (PT|) has two independent components: 

^ = (26) 
a} + lad 

= '^P' (27) 

where the generalized energy density 9" in this approximation is equal to the energy 
density e and the generahzed pressure is given by 

„2 



a 

The dot denotes the difierentiation with respect to time. 
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From the conservation law for the generahzed quantities: 

ts + 3H{Q + p) = 0, (29) 
where if = ^ is the Hubble parameter, we obtain the scaling law for the energy density: 

where Eq is a positive constant proportional to the total energy of matter within a sphere 
of radius a in the limit — ^ 0. In the course of time, a increases but the right-hand 
side of flHUj) never becomes negative since the assumption nT 3> a ceases to hold when 
_ ^ a3/2(t), we obtain 

^(6r = 3«:X^-aV. (31) 

For a very late universe kT <^ a, and the Einstein equations in the second 
approximation of a small quantity ^ are |2Qj 

- 2^9t^u = ■^K,T^u{^ + + ^9t,u, (32) 

where A = || — The coupling between matter and the gravitational field is 

decreased by a factor which tends to 3/4 as T — >■ 0. Equation (jH^ can be rewritten as 

2 ~ 3a 

From the same equation it also follows that the energy-momentum tensor in the Einstein 
frame is not covariantly conserved: || 

iiu 3 _|_ kT ' V"-*^/ 

4 4a 

unless T^y — \Tg^i, = 0.^ In the Jordan frame such a tensor is always conserved [2011221, 
although |TT] arrives at the different conclusion. 

The generalized energy density in the late FLRW universe with the metric ()24|) is 
given by 

^ 3e 3a 3Ke^ , , 

and the generahzed pressure is 

3a ne^ , 
IbK Iba 

From the conservation law (j29p [201 obtain the scaling law for the energy density: 

ea^ = Eooe"^, (38) 

II It might seem that this tensor is not conserved for a = 0. In this case, however, T — Q since kT <^ a. 
^ Covariant differentiation of the general equation gives 

which for an isotropic universe yields (|29|) . 
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where E^o is a positive constant. In the hmit of very small e, the above relation simplifies 
to the usual expression for nonrelativistic matter, ea^ = E^, where Eoo has the meaning 
of a quantity proportional to the total energy of matter within a sphere of radius a. 
Similarly, we obtain the energy density dependence on time: 

In^ - iMi^ = (39) 
eo 6a 4 

where time is measured from the instant at which e is equal to a given value eo- After 

some mathematical manipulations, we arrive at the following expression for the scaling 

factor as a function of time in the quadratic approximation: 

a3(t) = :^(e3^^c.t _ ^ + ^I^f^e-'""-'), (40) 
eo a 8a^ 



where H^a = c-y/a/16 [201 • The late universe approaches a de Sitter spacetime 
exponentially fast, as in the Jordan frame [H]. 

We note that the energy density of matter scales like (nonrelativistic matter) in 
both considered limits, Ke ^ a and ne ^ a. Therefore, there exists the largest deviation 
from such a scaling in the present epoch when and a are on the same order. 

4. The present universe in R+ 1/R gravity 

From the 00 component of equation we can obtain the Hubble parameter for a fiat 
FLRW universe filled with dust, as a function of c 



For the Lagrangian (fT^ . we find 



Hi<l>) = ^ J- - K (42) 
Similarly, we derive the deceleration parameter 5' = — ft using the 11 component: 

For the Lagrangian (fT^ . we obtain 

The time dependence of q is determined by combining equation ()44|1 with the function 
= 0(t) which is an integral of |20j : 

(45) 



We also find the expression for the apparent matter density parameter Q 



, 0oLX0o)-2L(0o) 

VIm = 2L (0o) (46) 

00^^ (0o) - 3L(0o) 
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where ec = is the critical energy density and the subscript refers to the present 
time. For the Lagrangian (fT^ . this parameter becomes 

Substituting 00 = 0o(q'o) from equation pij) to (jlTj) gives 

= (48) 
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or 



go 2 



^ (Qm - ^nl, - 5fiM + 4) . (49) 



For the early universe |0| ^ a and q ^ ^. For the late universe (f) —a and 
g — s> — 1. The transition from the decelerated (g > 0) to accelerated (g < 0) phase 
occurs at g = 0, which takes place when equals 



5 




(50) 

The Hubble parameter and the energy density of matter at the time of the deceleration- 
to-acceleration transition are given, respectively, by 

Ht = J^a, (51) 




We see that the energy density at this time is on the order of ^, i.e. it scales with the 
largest deviation from nonrelativistic dust. 

In the Jordan frame, the energy density scales like a~'^ which is equivalent to (l+z)^, 
where z is the redshift. The energy density in the Einstein frame is related to that in 
the Jordan frame by the conformal factor L'{(j)). Therefore, 

^ = ^°(1 + ^)'^- (53) 

^ (0o) 

Substituting here e(0) from equation (fT^ gives 
0L-(0) - 2L(0) [L'{<Po)f ^ 
0o^'(0o) - 2L(0o) [L'i<l>W ^ ^ ' ^ ' 

from which we obtain = (j){z, 0o)- In the case of the Lagrangian (jl4j) . is a solution 
of a quintic equation. From the redshift dependence of we can derive the relations 
H = H{z) and g = q{z) using equations and (jSj), respectively. 



The observed value of the present matter density parameter is VLm = 0.29± 



0.05 
0.03 



Substituting this value into equation gives the prediction for the apparent 

deceleration parameter: 

go = -0.67±°:°^, (55) 
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which agrees with the observed go = — 0.74±0.18 j21]-''' The deceleration-to-acceleration 
transition happened in the past and the universe expands with a positive acceleration. 
The corresponding value of 0o is determined by equation ()44|1 : 

00 = (-1.05±0.01)a. (56) 

From equation (@2I), we find 

a = (12.5±rr)^. (57) 

The observed value Hq = 71 ± A^^^ ^ leads to 

a = {7.35±lfj) X IQ-^^m-^. (58) 

The deceleration-to-acceleration redshift Zt can be obtained from equations (jl9p . 
(jSO)), (fS^ and (jSni)- Substituting the numerical value for 0o from yields 

Zt = 0.86 ±°:?o • (59) 

This value overlaps with the ACDM value of zt = 0.71±o;^g* but disagrees with the 
observed Zt = 0.46 ± 0.13 [5^. The extraction of Zt from the data appears, however, 
to be much less robust than the extraction of go |2Ej- Therefore, we need stronger 
measurements of Zt to verify if the R+l/R gravity is a viable explanation of current 
cosmic acceleration. 

The last non-dimensional parameter in our study is qi, the coefficient in the Taylor 
expansion of q{z) around z = 0: q{z) = qo + qiz. We use qi = §f U=o = to 
obtain 

^ L'{^L''~^LL"-LL'){^L'-2L) 

(0L'-3L)2(4LL"-L'2 -0L'L") *=<>o ^ ' 

from equations (|43p and (j54p . The Lagrangian (jl4p and the current value of (j56p give 
gi = 0.81±|ii^, (61) 
which slightly overlaps with the observed gi = 1.59 ± 0.63 

5. Summary 

Current cosmic acceleration can be explained by adding a. 1/R term to the Einstein- 
Hilbert Lagrangian. We used the metric-affine variational formalism and chose the 
Einstein conformal frame to be physical. In this frame we derived the field equations for 
the early and late matter-dominated FLRW universe, in the quadratic approximation 
with respect to small quantities. We showed that the largest deviation of the energy 
density e scaling from the usual scaling of nonrelativistic dust occurs in the present 
epoch. We did not give the exact form of the field equations for any value of e. 

The predictions for non-dimensional cosmological parameters in the i?+ 1/R gravity are independent 
of the value of a. 



For a flat universe in the ACDM model, zt is a simple function of flM- zt ~ 



2(1-Om) 



1/3 



1. 



The present universe in the Einstein frame, metric-affine R+l/R gravity 



10 



although it can be done analytically for the R + l/R Lagrangian. Instead, we derived 
the expressions for non-dimensional cosmological parameters to verify if the R+l/R 
gravity is compatible with observations. 

We found that the current deceleration parameter go derived from the apparent 
matter density parameter Qm is consistent with observations. There is also a tiny 
overlap between the predicted and observed values for the redshift derivative qi of 
the deceleration parameter. The predicted redshift of the deceleration-to-acceleration 
transition Zt agrees with that in the AC DM model but it is larger than the value 
estimated from SNIa observations. Since the robustness of the Zf measurements is 
weaker than that of qo, the question on the viability of the R+l/R gravity in the 
Einstein frame remains open. 
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